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Abstract 

The determination of the densest packings of regular tetrahedra (one of the five Platonic solids) 
is attracting great attention as evidenced by the rapid pace at which packing records are being 
broken and the fascinating packing structures that have emerged. We have discovered the densest 
known packings of regular tetrahedra with a density cj) = ]^\g = 0.855506 .... These packings are 
special cases of an analytical two-parameter family of dense periodic packings with four particles 
per fundamental cell that we have constructed here. From this family, we can recover a set of recent 
packing arrangements due to Kallus et al. with density <p = = 0.854700. . ., which has higher 
symmetry than our densest packings. We also describe a procedure that could lead to rigorous 
upper bounds on the maximal density of tetrahedron packings, which could aid in assessing the 
packing efficiency of candidate dense packings. 

PACS numbers: 61.50.Ah, 05.20.Jj 
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I. INTRODUCTION 



Dense packings of nonoverlapping solid objects or particles are ubiquitous in synthetic 
and natural situations. Packing problems arise in technological contexts, such as the pack- 
aging industries, agriculture (e.g., grains in silos), and solid- rocket propellants, and underlie 
the structure of a multitude of biological systems (e.g., tissue structure, cell membranes, 
and phyllotaxis). Dense particle packings are intimately related to the structure of low- 
temperature states of condensed matter, such as liquids, glasses and crystals [l-4|. In the 
last decade, scientific attention has broadened from the study of dense packings of spheres 
(the simplest shape that does not tile Euclidean space) {5 -14] to dense packings of nonspher- 
I5M17I] . A basic characteristic of a packing in d- dimensional Euclidean space 



ical particles 



M. d is its density <p, defined to be the fraction of M. d that is covered by the particles. 

A problem that has been of great scientific interest for centuries is the determination of 
the densest arrangement (s) of particles that do not tile space and the associated maximal 
density 4> m ax- For generally shaped particles, finding the densest packings is notoriously 
difficult. This salient point is summarized well by Henry Cohn [3] who recently remarked, 
"For most grain shapes we cannot guess or even closely approximate the answer, let alone 
prove it, and it is difficult to develop even a qualitative understanding of the effects of grain 
shape on packing density." Until recently, very little was known about the densest packings 
of polyhedral particles. The difficulty in obtaining dense packings of polyhedra is related to 
their complex rotational degrees of freedom and to the non-smooth nature of their shapes 



19 



20]. 



Recently, we set out to attempt to determine the densest known packings of the Platonic 
and Archimedean solids {19I . 20]. It was shown that the central symmetry of the majority 
of the Platonic and Archimedean solids distinguish their dense packing arrangements from 
those of the non-centrally symmetric ones in a fundamental way. (A particle is centrally 
symmetric if it has a center C that bisects every chord through C connecting any two bound- 
ary points of the particle; i.e., the center is a point of inversion symmetry.) The tetrahedron 
is the only Platonic solid that lacks central symmetry, an attribute that geometrically frus- 
trates it to a greater degree than the majority of the remaining solids in this set that do 



not tile space 



20J. A number of organizing principles emerged in the form of conjectures for 



polyhedra as well as other nonspherical shapes. In the case of polyhedra, the following two 
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are particularly relevant: 



• Conjecture 1: The densest packings of the centrally symmetric Platonic and 
Archimedean solids are given by their corresponding optimal Bravais lattice packings. 

• Conjecture 2: The densest packing of any convex, congruent polyhedron without 
central symmetry generally is not a (Bravais) lattice packing, i.e., set of such polyhedra 
whose optimal packing is not a Bravais lattice is overwhelmingly larger than the set 
whose optimal packing is a Bravais lattice. 

Conjecture 1 is the analog of Kepler's sphere conjecture for the centrally symmetric Pla- 
tonic and Archimedean solids. In this sense, such solids behave similarly to spheres in that 
their densest packings are lattice arrangements and (except for the cube) are geometrically 



frustrated like spheres. Conjecture 2 has been shown by Conway and Torquato [21] to be true 
for both the tetrahedron and Archimedean non-centrally symmetric truncated tetrahedron, 
the latter of which can be arranged in a "uniform" non-Bravais lattice packing with density 
at least as high as 23/24 = 0.958333 .... (A uniform packing, defined more precisely below, 
has a symmetry that takes one tetrahedron to another.) Indeed, it was this investigation 
that has spurred the flurry of activity in the last several years to find the densest packings of 
tetrahedra. There have been many twists and unexpected turns since 2006 that have led to 
the densest known packings of tetrahedra that we report here with = ^|||| = 0.855506 . . ., 
which is the highest density obtained from our two-parameter family of constructions de- 
scribed below. Therefore, to place our present results in their proper context, it is instructive 
to review briefly the developments since 2006. 

First, we note that the densest Bravais-lattice packing of tetrahedra (requiring one tetra- 
hedron per fundamental cell such that each tetrahedron in the packing has the same orienta- 
tion as the others) has a packing fraction <j> = 18/49 = 0.367 . . . and each tetrahedron touches 
14 others 25j. Conway and Torquato [21] showed that the densest packings of tetrahedra 
cannot be Bravais lattices by analytically constructing several such packings with densities 
that are substantially larger than 0.367. (A non-Bravais lattice packing contains multiple 
particles, with generally different orientations, per fundamental cell, which is periodically 
replicated in M. d .) One such packing is a "uniform" packing with density = 2/3 and two 
particles per fundamental cell. The so-called "Welsh" packing has a density = 0.708333 . . . 
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TABLE I: A brief summary of the dense non-lattice packings of tetrahedra. The name of the 
packing is given along with the year that it was discovered. Here <ft is the packing density and N 
is the number of tetrahedra per fundamental cell. 



Packing 


Year 


<t> 


N 


Uniform I [21] 


2006 


| = 0.666666 . . . 


2 


Welsh [21] 


2006 


±f = 0.708333 . . . 


34 


Icosahedral [21] 


2006 


0.716559... 


20 


Uniform II [22] 


2009 


139+ 3 ^° = 0.719488.. 


2 


Wagon Wheels [23] 


2008 


0.778615... 


18 


Improved Wagon Wheels [19] 


2009 


0.782021 . . . 


72 


Disordered Wagon Wheels [20] 


2009 


0.822637... 


314 


Ring Stacks [24] 


2009 


0.8503... 


82 


Uniform III [22] 


2009 


= 0.854700 . . . 


4 


Dimer-Uniform 


2009 


14319 = 0.85 5 5 06 . . . 
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and 34 particles per fundamental cell. Yet another non-Bravais lattice packing with den- 
sity = 0.716559... is based on the filling of "imaginary" icosahedra with the densest 
arrangement of 20 tetrahedra and then arranging the imaginary icosahedra in their densest 
lattice packing configuration. The densities of both the Welsh and Icosahedral packings 
can be further improved b y c ertain particle displacements 21] . Using imperfect "tetrahe- 
dral" dice, Chaikin et al. 26j] experimentally generated jammed disordered packings with 



~ 0.75. Employing physical models and a computer algebra system, Chen |23[] discovered 
a remarkably dense periodic arrangement of tetrahedra with = 0.7786 . . ., which exceeds 
the density (0 ma z = 7t/a/18 = 0.7404. . .) of the densest sphere packing by an appreciable 
amount. We have called this the "wagon-wheels" packing 19j, |20|. 

Torquato and Jiao [19J devised and applied an optimization scheme, called the adaptive- 
shrinking-cell (ASC) method, that used an initial configuration based on the wagon-wheels 
packing to yield a non-Bravais lattice packing consisting of 72 tetrahedra per fundamental 



cell with a density = 0.782. . . 19j. Using 314 particles per fundamental cell and start- 



ing from an "equilibrated" low- density liquid configuration, the same authors were able to 



improve the density to = 0.823 . . . 20j. This packing arrangement interestingly lacks long- 
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range order. Haji-Akbari et al. [24J numerically constructed a periodic packing of tetrahedra 
made of parallel stacks of "rings" around "pentagonal" dipyramids consisting of 82 particles 



per fundamental cell and a density = 0.8503 .... More recently, Kallus et al. [22j found 
a remarkably simple uniform packing of tetrahedra with high symmetry consisting of only 
four particles per fundamental cell with density <p = ^ = 0.854700 .... Table E summarizes 
some of the packing characteristics of the non-Bravais lattice packings of tetrahedra. 



II. TWO-PARAMETER FAMILY OF DENSE PACKINGS OF TETRAHEDRA 

Inspired by the work of Kallus et al. j^j], we have applied the adaptive-shrinking-cell 
(ASC) optimization scheme to examine comprehensively packings with a considerably small 
number of particles per fundamental cell (from 2 to 32) than we have used in the past. The 
ASC scheme employs both a sequential search of the configurational space of the particles 
and the space of lattices via an adaptive fundamental cell that deforms and shrinks on 
average to obtain dense packings A dense packing with 8-particle basis that emerged from 
this numerical investigation suggested that it was composed of two very similar fundamental 
cells, each containing 4 particles. Using one of the 4-particle basis configurations, we were 
able to find packings with density = 0.8551034... that exceeded the highest density 
packings with <fi = 100/117 = 0.854708... constructed by Kallus et al. Even though our 
packings possess a type of point inversion symmetry, they are not as symmetric as the densest 



packings reported in Ref. [22j], as we now explain. The four tetrahedra in the fundamental 
cell in our dense numerically generated packings formed two contacting "dimers" . A dimer is 
composed of a pair of regular tetrahedra with unit edge length that exactly share a common 
face. The compound object consisting of the two contacting dimers possesses point inversion 
symmetry, with the inversion center at the centroid of the contacting region on the faces. A 
Bravais lattice possesses point inversion symmetry about the lattice points and the centroids 
of the fundamental cells. By placing the symmetry center of the two-dimer compound on 
the centroids (or the lattice points), we construct packings that generally possess point 
inversion symmetry only about the symmetry centers of the two-dimer compound. We call 
such structures dimer-uniform packings, since the inversion symmetry acts to take any dimer 
into another. Such packings should be distinguished from the more symmetric uniform 
(or transitive) packings of tetrahedra in which the symmetry operation acts to take any 
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tetrahedron into another, such as the ones found in Refs. [2l| and 22] (see Tabled]). The 
latter have almost as much symmetry as a Bravais lattice, except that the centroids of the 
particles are not just characterized by simple translational symmetry. 

We then set out to obtain analytical constructions based on our numerical packings that 
relax the symmetry conditions on the contacting dimers. In particular, we orient the 3-fold 
rotational symmetry axis of one of the dimers in an arbitrary direction (say the ^-direction of 
a Cartesian coordinate system), and then fix the origin of the lattice vectors at the centroid 
of this dimer. (The centroid is located at the center of the contacting faces of the two 
tetrahedra that comprise the dimer.) Then we place the second dimer in contact with the 
first one such that there is a center of inversion symmetry that takes one dimer to the other, 
which implies face-to-face contacts between the two dimers. 

The problem of determining the analytical constructions then amounts to determining 
12 equations for the 12 unknowns. Nine of the 12 unknowns arise from the three unknown 
lattice vectors, each of which contains three unknown components. The other 3 unknowns 
derive from the components of the centroid of the second dimer. 

In particular, we let the centroid of the dimer at the origin denoted by v G , and the centroid 
of the other dimer denoted by v c . The vertices of the two dimers at v c and v c are given by 

r ^ = (i 57s' °)> r * = H> 57s' °)> r ° = (°> -73' °)> TD = (°> °> yi } ' r£ = (0 ' °' " and 
r* A = -r A + v c , r* B = -r B + v c , v* c = -r c + v c , r* D = -r D + v c , = -i£; + v c , respectively. 

In addition, let the lattice vectors be v 1; v 2 , and v 3 . The 12 components of the four vectors 

(vj, i=c, 1,2,3) are the aforementioned unknowns that are related to each other through the 

nonoverlapping conditions. 

In our packings, each dimer has 8 face-to-face contacts and at least 2 edge-to-edge con- 
tacts. Among these 20 contacts of the two dimers in the fundamental cell, there are 8 
independent face-to-face contacts and 1 independent edge-to-edge contact, which reduces 
the number of variables for the packing from 12 to 3. 

A face-to-face contact requires that the projection of the vector distance between the 
centroids of the two dimers on the contacting face normal is a constant. The 8 independent 
face-to-face contacts are between the dimer pairs with the centroids at {v c , v Q }, {v c , v Q +vi}, 
{v c , v + v 2 }, {v + v 3 , v c }, {v + v 3 , v c - v 2 }, {v + v 3 , v c + vi - v 2 }, {v , v c - Vi + v 3 }, 
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and {v D , v c + vi — V2 — 2V3}. The contact between dimer pairs at {vj, Vj} requires 

(V< - Vj-) • Ilij = — , (1) 

where is unit outward contacting face normal of the dimer at Vj. 

The edge-to-edge contact requires that the projection of the vector connecting the corre- 
sponding ends of two edges on the common perpendicular line of the two edges equals zero, 



[T A - (T B - V! + v 2 + v 3 )] • 1 Q = 0, (2) 

where 1 = (r A - r D ) x (r B - r D ). 

Moreover, there are two independent nonoverlapping conditions obtained from 4 possible 
edge-to-edge contacts between neighboring particles, i.e., 

[r c - (r B - V! + v 3 )] • li > 0, (3) 

[r c ~ (t d + v 3 )] • h > 0, (4) 

where li = (r B - r D ) x (r E - r c ) and 1 2 = (r D - r A ) x (r c - r E ). 

Furthermore, there are two additional nonoverlapping conditions given by 4 potential 
vertex-to-face contacts, i.e., 

v 2 • ni = a/!, (5) 
where rix = (^, — — |) is unit outward normal of the contacting face, 

(vi - v 2 ) ■ n 2 = y|, (6) 

and n 2 = (— 3 , — , |) is unit outward normal of the contacting face. 

The edge-to-edge and vertex-to-face contacts are realized when the equality holds in the 
above conditions (J2J)-(l6]) . However, these contacts can not be realized simultaneously in 
general. In particular, Eqs. (J2K, ©, (ED) and Eqs. (CED, (PJ), @, (J5J) provide two sets of 
equations that lead to a family of dense tetrahedral packings whose structures are determined 
by two parameters a and b, with the corresponding density dependent only on the parameter 
a. In other words, each density is associated with a spectrum of different packing structures. 

Finally, we arrive at an explicit expression for the packing density: 



221 ] are recovered 



FIG. 1: (color online). The density <f) surface of our family of tetrahedral packings as a function of 
the two parameters a and b. As explained in the text, the thick red lines show two sets of tetrahedral 
packings with distinct structures but with the same density. The two red points correspond to the 
densest known two tetrahedral packings. The packings found by Kallus et al. 
from our two-parameter family (thick blue line). 

where Vr = V2/12 is the volume of a regular tetrahedron with unit edge length. Substituting 
the lattice vectors expressed in terms of the two parameters a and b into (j7|) yields 

100 

" 117 - 240a" (8) 

where a G ( — jfo ? ifo)- ^ ^ s important to note that for each a ^ 0, there are two sets of 
packings of tetrahedra, each with distinct structures but possessing the same density (as 
shown in Fig. [I] by the thick red lines). For — j|g < a < 0, the packings are specified by 

A a 4 3^2 [2 . 

Vc = ( 5 + 3'-57!'~57! + V3 a) ' 

vi = (-a, — — , y -a), 

V2 = ^ + a-b,~^ + V3a,-^=-^=a + V6b), 

,7 , V3 r 9 10 

v 3 = (---2a + b,---V3a,-— + -a-V6b), 



(9) 
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(a) 



(b) 



FIG. 2: (color online). Two different configurations of the densest known packings of four tetrahe- 
dra (two dimers) within their corresponding rhombohedral fundamental cells: (a) Here a = —3/140 
and 6 = 0, and the fundamental cell is colored yellow and its boundaries are colored red. (b) Here 
a = 3/140 and 6 = 0, and fundamental cell is colored blue and its boundaries are colored red. The 
two packings shown in (a) and (b) are only slightly different from one another. Specifically, the 
difference between the coordinates of the centroids that are not at the origin is (1/70, 0, \/6/70). 
The differences between the shapes of the rhombohedral fundamental cells in the two cases are 
readily apparent. 

where < b < g ° a - For < a < the packings are specified by 

,1 a 4 3^2 [2 , 

Vc = ( H + 3'-o7!'"o7! + V3 a) ' 

vi = {-a, — — , y -a) 

/3 „ , V3 at 3 10 /-. 
v 2 = (- - 2a - b, -— + VSa, -— + —a + V66), 

7 , f- 9 



(10) 



V3 = (-- + a + 6,---V3a,-— --a-V66), 
where < b < 

The densest packings are associated with a = — b = and a = 6 = 0, possessing 
a density 4> ma x = jf|f§ = 0.855506.... In each set, there is a unique packing structure asso- 
ciated with (fimax (shown as the red points in Fig. CD), instead of a spectrum of structures. 
Two different configurations of the densest known packings of four tetrahedra within their 
corresponding rhombohedral fundamental cells are shown in Fug. [2j Figure [3] depicts peri- 





(a) (b) 

FIG. 3: (color online). This figure shows periodic replicates the densest known tetrahedral packings 
corresponding to Fig. [2] with 8 fundamental cells (2 along each lattice vector. The four tetrahedra 
within a fundamental cell are shown in blue: (a) a = —3/140 and 6 = 0. (b) a = 3/140 and 6 = 0. 
Note that in (a), the dimer with centroid at vi is slightly shifted to the right with respect to the 
dimer at the origin; and in (b), the dimer with the centroid at vi is slightly shifted to the left with 
respect to the dimer at the origin. 

odic replicates of our optimal tetrahedral packings corresponding to those shown in Fig. [2] 
with 8 fundamental cells (2 along each lattice vector). 

At a = 0, there is only one set of packings with the same <fi = 100/117 = 0.85470..., 
whose structures are dependent on b (shown as the blue line in Fig. [1]). These packings 
reduce exactly to those discovered by Kallus et al., which possess relatively high symmetry. 
In particular, a = allows the centroids of the dimers, which are related to each other 
by an integer multiple of vi, to be perfectly aligned on two of the mirror image planes of 
the dimers simultaneously, which leads to additional two-fold rotational symmetry of the 
packing. This additional rotational symmetry, together with the point inversion symmetry, 
leads to uniform packings with respect to each tetrahedron (not just each dimer), i.e., the 
symmetry operation acts to take each tetrahedron into another. 
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III. TOWARDS UPPER BOUNDS ON THE MAXIMAL DENSITY 



The problem of determining upper bounds on the maximal density of packings of non- 
spherical particles is highly nontrivial, and yet such estimates would be indispensable in as- 
sessing the packing efficiency of a candidate dense packing, especially if tight upper bounds 
could be constructed. It has recently been shown that (f) max of a packing of congruent 
nonspherical particles of volume vp in M 3 is bounded from above according to 



< min 



, (ID 

v s Vl8 J 

3e inscribed in the nonspherical particle 
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20J. The upper bound (ITT]) will be 



where vs is the volume of the largest sphere that can 
and 7t/a/I8 is the maximal sphere-packing density 
relatively tight for packings of nonspherical particles provided that the asphericity 7 (equal 
to the ratio of the circumradius to the inradius) of the particle is not large. However, for 
tetrahedra, the asphericity is too large for the upper bound ( TTTT) to yield a result that is less 
than unity. 

One possible approach to obtaining nontrivial upper bounds is to attempt to generalize 



the idea that Rogers used to prove upper bounds on <p max for sphere packings [27|. The 
key concept is to consider a locally dense cluster of 4 contacting spheres in a tetrahedral 
arrangement and then prove that the fraction of space covered by the spheres within the 
tetrahedron joining the sphere centers is an upper bound on <p m ax- This can be done because 
one can triangulate any sphere packing to decompose it into generally irregular tetrahedra 
with vertices at sphere centers. The fact that the regular tetrahedron has the best density 
for any tetrahedron, then yields an upper bound for the density of any sphere packing. In 
the case of the non-tiling Platonic and Archimedean solids, a natural choice for the enclosing 
region associated with the cluster is its convex hull. 

For tetrahedra, we must identify the least densest local cluster with density that exceeds 
<t>max- A trivial choice is a dimer because the fraction of space covered by the dimer within 
its convex hull is unity. A nontrivial choice is a 5-particle "wagon-wheel" cluster as shown 
in Fig. HI The convex hull of the local packing of five tetrahedra sharing a common edge 
(a "wagon wheel" cluster) can be decomposed into five regular tetrahedra and five thin 
irregular tetrahedra. We assume the gaps between the regular tetrahedra are equal, i.e, the 
thin irregular tetrahedra are congruent. Note the regular tetrahedron shares two faces with 
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FIG. 4: (color online). The convex hull of five regular tetrahedra in a "wagon wheel" arrangement. 
The convex hull can be decomposed into five regular tetrahedra (shown in blue) and five thin 
irregular tetrahedra (shown in red). 

its neighboring irregular tetrahedra. Thus, the volume ratio is equal to the ratio of the 
corresponding heights of the regular and irregular tetrahedron, i.e., 

V T 1 
7 ~ ^(3cos 2 ff -1)' (12) 

where Vr and Vr„ is the volume of the regular and irregular tetrahedron, respectively. Thus, 
the density of this local packing, denned as the fraction space covered by the regular tetra- 
hedra within the convex hull is given by 

V T 1 



>w 



V T + V n 3yJ cos 2^ + (1 _y^N 



0.974857.... (13) 



Because this cluster is highly anisotropic and is an effectively "flat" object, it is reasonable 
to assume that it is not the least densest local cluster and therefore its local density is a 
gross overestimate of 4> ma x- 

Since tetrahedra are fully three-dimensional objects, it seems reasonable to assume that 
the least densest local cluster should be more "isotropic" than the highly anisotropic and 
effectively "two-dimensional" wagon-wheel cluster. One plausible choice for such a cluster 
consists of 20 tetrahedra sharing a common vertex (i.e., an icosahedral-like cluster). The 
convex hull of this cluster can be decomposed into 20 regular tetrahedra, 12 pyramids with 
pentagonal bases and 30 pyramid with rectangular bases (see Fig. [5]). We assume the gaps 
between tetrahedra are equal and thus the two types of pyramids are congruent. The volume 
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FIG. 5: (color online). The convex hull of 20 regular tetrahedra in an "icosahedral" arrangement. 
The convex hull can be decomposed into 20 regular tetrahedra (shown in blue), 12 pyramids with 
pentagonal bases (shown in yellow) and 30 pyramid with rectangular bases (shown in red). 

of the regular tetrahedron is Vp = and the volume Vp of the pyramid with pentagonal 



base is given by 



5 ,3tt. r , 
V P = — tan( — )L 2 
P 12 v 10 ; 



L 2 



4 cos 2 ff ' 



(14) 



and the volume Vr of the pyramid with rectangular base is given by 



V R = ^Ly[Y^L\ 
o 

where 



(15) 



and r 



L 



.2^2 



- 1 




(16) 



V / 5)/2 is the golden ratio. Thus, the local packing density is given by 



20V r 



- 2 v T + i2v P + s w- R = 0Mm55 - (17 > 

Note in the above calculations, we have assumed that the gaps between the tetrahedra 
are equal, which is sufficient to provide an estimate of the fraction of space covered by the 
cluster within its convex hull. However, the idea of Rogers to prove an upper bound for 
sphere packings cannot be used here because there is no analogous decomposition of space 
into irregular convex hulls shown in Figs. |4] or [5j A completely new idea is needed to prove 
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that the aforementioned estimates are bounds, but it seems plausible that they are correct 
bounds. If one could prove that these clusters are indeed locally denser than the globally 
densest packings, then the aforementioned estimates provide upper bounds on 4> max , but 
they cannot be sharp, i.e., they are not achievable by any packings. It is noteworthy the 
density <p = 0.855506 ... of our our densest packings is relatively close to this putative upper 
bound density estimate of 0.880755 .... 



IV. DISCUSSION 



For all of the small systems (including 2 to 32 particles) that we investigated using our 
numerical ASC scheme, the densest packings that emerged had a 4-particle basis. Our 
analytical constructions indicate that for such small systems, the highest density packings 
we found could be optimal. However, there is no reason to believe that denser packings 
could not be discovered by carrying out exhaustive searches to determine the globally maxi- 
mal densities of packings with successively larger numbers of particles per fundamental cell. 
Previous numerical studies have indicated that dense packings may have a large number of 
particles per fundamental cell arranged in a complex fashion, e.g., the "disordered wagon- 
wheels" packing with = 0.822637 . . . [2 Of and the "ring stacks" packing with = 0.8503 . . . 



241 ] . However, it is plausible that such packings are in fact only locally optimal solutions 



and hence the numerical techniques used to obtain them are incapable of extricating them- 
selves from these "trapped" regions of configuration space to find denser and more ordered 
structures due to the intrinsic geometrical frustration of the tetrahedron. This may also 
call into question claims made by Haji-Akbari et al. 24| that their packings, which are 
characterized by an effective "quasicrystal-like" plane, are true thermodynamic equilibrium 
phases of tetrahedra, especially at high densities. Instead, our highest-density constructions 
suggest that uniform periodic packings with a 4-particle basis (or even some yet unknown 
denser periodic packing) and unjammmed, lower-density countyerparts could be the stable 
phases at such high densities. If the latter is correct and the putative "quasicrystal-like" 
phase truly exists at intermediate densities, then it is hard to imagine how a quasicrystal-like 
structure of tetrahedra under quasi-static compression (densification) could rearrange to a 
structure with a more ordered periodic arrangement with higher symmetry. However, it is 
difficult to draw such definitive conclusions without further study. 
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Although there could still be tetrahedral packings denser than our constructions, it ap- 
pears that all of the evidence thus far points to the fact that the densest tetrahedral pack- 

1J2I, 23 



ings cannot possess very high symmetry 



24j due to the lack of central symmetry 



of a tetrahedron and because tetrahedra cannot tile space [20|]. Indeed, our densest uni- 
form 4-particle-basis packings found in the present paper improved upon the best analogous 
packings of Kallus et al. 22j by relaxing the two-fold rotational symmetry constraints they 
imposed. 



Acknowledgments 

We are grateful to Henry Colin, Yoav Kallus and John Conway for helpful discussions. 
S. T. thanks the Institute for Advanced Study for its hospitality during his stay there. This 
work was supported by the Division of Mathematical Sciences at the National Science Foun- 
dation under Award Number DMS-0804431 and by the MRSEC Program of the National 
Science Foundation under Award Number DMR-0820341. 



[1] J. D. Weeks, D. Chandler, and H. C. Andersen, J. Chem. Phys. 54, 5237 (1971). 
[2] R. Zallen, The Physics of Amorphous Solids (Wiley, New York, 1983). 

[3] P. M. Chaikin and T. C. Lubensky, Principles of Condensed Matter Physics (Cambridge 

University Press, New York, 1995). 
[4] S. Torquato, Random Heterogeneous Materials: Micro structure and Macroscopic Properties 

(Springer- Verlag, New York, 2002). 
[5] C. F. Gauss, Gottingsche Gelehrte Anzeigen (1831), see also J. reine angew. Math., vol. 20, 

1840, 312-320. 
[6] J. D. Bernal, Nature 185, 68 (1960). 

[7] J. H. Conway and N. J. A. Sloane, Sphere Packings, Lattices and Groups (Springer- Verlag, 

New York, 1998). 
[8] A. J. Liu and S. R. Nagel, Nature 396, 21 (1998). 

[9] S. Torquato, T. M. Truskett, and P. G. Debenedetti, Phys. Rev. Lett. 84, 2064 (2000). 
[10] C. S. O'Hern, S. A. Langer, A. J. Liu, and S. R. Nagel, Phys. Rev. Lett. 88, 075507 (2002). 



15 



[11] H. Cohn and N. Elkies, Annals Math. 157, 689 (2003). 

[12] T. C. Hales, Ann. Math. 162, 1065 (2005). 

[13] S. Torquato and F. H. Stillinger, Experimental Math. 15, 307 (2006). 

[14] H. Cohn and A. Kumar, Ann. Math. 170, 1003 (2009). 

[15] U. Betke and M. Henk, Comput. Geom. 16, 157 (2000). 

[16] A. Donev, F. H. Stillinger, P. M. Chaikin, and S. Torquato, Phys. Rev. Lett. 92, 255506 
(2004). 

[17] Y. Jiao, F. H. Stillinger, and S. Torquato, Phys. Rev. E 79, 041309 (2009). 

[18] H. Cohn, Nature 460, 801 (2009). 

[19] S. Torquato and Y. Jiao, Nature 460, 876 (2009). 

[20] S. Torquato and Y. Jiao, Phys. Rev. E 80, 041104 (2009). 

[21] J. H. Conway and S. Torquato, Proc. Nat. Acad. Soc. 103, 10612 (2006). 

[22] Y. Kallus, V. Elser, and S. Gravel (2009), arXiv:0910.5226. 

[23] E. R. Chen, Discrete Comput. Geom. 40, 214 (2008). 

[24] A. Haji-Akbari, M. Engel, A. S. Keys, X. Zheng, R. G. Petschek, P. Palffy-Muhoray, and S. C. 

Glotzer, Nature 462, 773 (2009). 

[25] D. J. Hoylman, Bull. Amer. Math. Soc. 76, 135 (1970). 

[26] P. M. Chaikin, S. Wang, and A. Jaoshvili, in Am. Phys. Soc. March Meeting (2007). 

[27] C. A. Rogers, Proc. Lond. Math. Soc. 8, 609 (1958). 



16 



